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1 Introduction 



It has now been firmly established that neutrinos have tiny masses and that oscillation occurs 
between different flavors. The standard paradigm that explains the small masses is the seesaw 
mechanism [1], which generates an effective dimension-5 operator 0\ = (LLHH)/M, suppressed 
by the mass scale M of the heavy right-handed neutrino. Oscillation data suggests that in this 
scenario M ~ 10 14 , which is well beyond the reach of foreseeable experiments for direct scrutiny. 
One has to rely on other indirect hints, such as charged lepton flavor violation in a supersymmetric 
context, in order to falsify this theory. 

An interesting alternative to the high scale seesaw mechanism is to induce small neutrino 
masses at the loop level. The smallness of neutrino masses can be understood as originating 
from loop and chirality suppression factors. The simplest among this class of models is the Zee 
model [3] where neutrino masses are induced as one-loop radiative corrections arising from the 
exchange of charged scalar bosons. The effective operator in this model is LLLe c H/M. To 
convert this to neutrino mass, a loop diagram is necessary. 1 In a second class of models, neutrino 
masses arise as two-loop radiative corrections via the exchange of singly and doubly charged 
scalars [5, 6]. The effective operator of these models is LLLe c Le c /M 2 . This model is compatible 
with neutrino oscillation data. Phenomenology of these models has been studied in Ref. [7]. 

A classification of low dimensional effective AL = 2 lepton number violating operators that 
can lead to neutrino masses has been given in Ref. [8]. Among these is an operator labeled O s : 

O s = L a ^WH p e aP (1) 

where a, (3 are SU (2) L indices, with the family indices suppressed. It is the purpose of this paper 
to develop a renormalizable model that generates 8 - We will show that such a model can be 
consistently constructed, with a variety of testable predictions. 

Operator Og is most directly induced by the exchange of scalar leptoquarks. 2 For neutrino 
mass generation the two quark fields and the charged lepton field in Og will have to be removed, 
which implies that the masses arise at the two-loop level. The order of magnitude of the induced 
neutrino masses is 

m t m h m T \iv 

where /i is a dimensionful coefficient of a cubic scalar coupling, and v = 174 GeV is the electroweak 
vacuum expectation value (VEV). In order to generate m v ~ 0.05 eV, it is clear that M L q must 
be of order TeV, which would be within reach of the LHC. 



1 This model is now excluded by neutrino oscillation data (see Ref. [4]). 

2 Neutrino mass generation at the one- loop level by leptoquark exchange has been studied in Ref. [9] . Neutrino 
mass generation at the two-loop level in the supersymmetric standard model with specific i?-parity violating 
couplings has been studied in Ref. [10]. 
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TeV scale leptoquarks can mediate a variety of flavor changing processes. We have examined 
these constraints and found the neutrino mass generation to be self-consistent. Our findings 
include several interesting features: (i) /j, — > may be suppressed, because of a GIM-like 
cancelation, but fi — > 3e and fi — e conversion nuclei are within reach of the next generation 
of experiments, (ii) There is a new CP violating contribution mediated by leptoquark box 
diagrams in B s —B s mixing, which can nicely fit the recently reported dimuon anomaly by the D0 
collaboration, (iii) The la discrepancy between theory and experiment in the leptonic decay of 
Df mesons can be explained in the model, owing to new contributions from the leptoquarks. (iv) 
Muon g — 2 receives new positive contributions which can resolve the theoretical anomaly there. 

(v) Neutrino oscillation parameters can be probed in the branching ratios of the leptoquarks. 

(vi) The leptoquark masses are constrained to be less than a few TeV, which should make them 
accessible to the LHC. 

This paper is organized as follows. In Section 2 we discuss the model leading to two-loop 
neutrino mass generation and obtain the constraints placed on the oscillation parameters. The 
constraints from rare processes such as n~ — > e _ 7 and fx~ — > e + e~e~ are presented in Section 3. 
In Sec. 4 we discuss the collider signals of leptoquarks, and in Sec. 5 we conclude. 

2 Model of two— loop neutrino mass generation 

In this section we present our model of two-loop neutrino mass generation, and derive restric- 
tions on the model parameters from neutrino oscillation data. Constraints from rare processes, 
discussed in Sec. 3, will be used to demonstrate the viability of the model and its predictions for 
neutrino oscillations. 



The gauge symmetry of our model is the same as the standard model (SM), SU(3) C x SU(2)l x 
U(1)y- In addition to the SM Higgs doublet H(l, 2, 1/2), the scalar sector consists of the following 
leptoquark multiplets: 



In general, addition of leptoquarks into the theory can cause baryon number (B) violating in- 
teractions, we forbid them by assuming that B is globally conserved. The leptoquarks have the 
following Yukawa interactions: 



2.1 Model 




X 



V3 -(3,1,-1/3). 



(3) 



Yukawa 



(4) 
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Here i, j — 1 — 3 are family indices and a, (3 are SU{2) L indices. Note that these Yukawa 
couplings conserve both baryon number and lepton number (L), as can be seen by assigning 
(B,L) charges of (1/3,-1) to f2 and (1/3,1) to x _1//3 - The couplings Y^u^d^x*, allowed by 
the gauge symmetry are forbidden by B, and the couplings F^LiQjx* , allowed by the gauge 
symmetry as well as B are forbidden by lepton number symmetry, which is assumed to be 
broken only by soft terms. 3 This breaking of L by two units occurs softly via a cubic term in the 
scalar potential: 

V = fj, ^H X ~ 1/3 + h.c. . (5) 

The simultaneous presence of Eqs. (4) and (5) would imply that neutrino masses will be generated 
at the loop level, as they lead to to the effective dimension 7 operator (Ld c )(u c e c )H [8], once the 
heavy leptoquark fields are integrated out. 

The Lagrangian relevant for neutrino mass generation in component form is 

C v = ^MK 1/3 - e ^ 2/3 ) + F ^XjX- 1/3 - ^- 2 % + + W % ) X - 1/3 + h.c. (6) 

Once the neutral component of the SM Higgs doublet acquires a vacuum expectation value (VEV) 
v = 174 GeV, the cubic term in the scalar potential will generate mixing between cu' 1 ^ 3 and x -1 ^ 3 
leptoquarks, with a mass matrix given by 

Mlo=( ml • (7) 

\ if v m l x J 

The parameter /x can be made real by redefining the leptoquark fields. We diagonalize this matrix 
to obtain the leptoquark mass eigenstates through 



x -i/s 

where cq = cos 9, sq = sin 9 with the angle 9 given by 




(8) 



2ut> . . 

tan2fl= - £ 9 



The squared masses of X l ^ 3 are given by 

^i,2 = l[™l + m 2 x T yJ(ml-m* x y + 4fiW ] (10) 



with Ml < M|. 



3 The couplings of course do not mediate rapid proton decay, however, their simultaneous presence with 
the couplings of Eq. (4) would lead to severe restrictions on F{j , since the successful V — A structure of the SM 
will then be drastically altered [11]. Although not essential, we prefer to set these F(j couplings to zero by L 
symmetry. 
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W H~ 

(a) W exchange (b) H exchange 



Figure 1: The two-loop diagrams contributing to neutrino mass generation. 

2.2 Two— loop neutrino mass diagrams 

Since lepton number violation occurs only when all three terms of Eq. (6) are simultaneously 
present, neutrino masses are generated in the model only at the two-loop level. The relevant 
diagrams, which involve the exchange of leptoquarks and a W ± boson, are shown in Fig. 1. We 
have evaluated these diagrams in the Feynman gauge. In this gauge, the unphysical charged 
Higgs boson exchange has to be kept. Interestingly, this set of charged Higgs diagrams add up 
to zero for the neutrino mass. The underlying reason for this result is that in the SM, both 
the up-type quark masses and the down-type quark masses are generated from the same Higgs 
doublet. Consequently the charged Higgs boson Yukawa couplings to the up-type quark has a 
relative minus sign compared to its Yukawa couplings to the down-type quarks, as given in the 
interaction Lagrangian 

C H ± = ~ [ur M u V d L H- -u L V M d d R H~ + h.c] (11) 

v 

Here V is the CKM matrix, and M u and M d are the diagonal mass matrices for the up and down 
quarks. When combined with Fig. lb, this relative minus sign implies zero net contribution to 
the neutrino mass from the charged Higgs boson exchange. In Fig. lb, the vertex u c l d c k H + is a 
sum of two contributions, one where chirality flip occurs in the u\ line, and the other where it 
occurs in the d\ line. These two contributions exactly cancel, see Eq. (11). 

Straightforward evaluation of the leptoquark-H^ exchange diagrams gives the neutrino mass 
matrix as: 

{M u )ij = m [MDMV^kiWtiF^D^ + (D £ ) i (F*) il (D u ) l V lk {D d ) k (Y T ) kj ] I jkl . (12) 

Here D u ^i are the (normalized) diagonal mass matrices for up quarks, down quarks, and charged 
leptons: 



J I I I I I I I I I I I I L_ 

0.2 0.4 0.6 0.8 



i I i i i I i i i I i i i I i i i I i i i I i i 

1.0 1.2 1.4 1.6 1.8 2.0 

M 1 (TeV) 



Figure 2: The exact loop integral with top quark inside the loop as a function of leptoquark masses Mi and 
M 2 . 



The overall scale m is given by 



m 



Cg 2 sin 26\ / m t mbm 7 
(16tt 2 ) 2 ) V Ml 



(14) 



where M\ is the lighter of the two charge —1/3 leptoquark mass and C = 3 is a color factor. 
The dimensionless two-loop integral Ijki depends on the two leptoquark masses, the W boson 
mass, m w , and the up-type quark masses m Ul , down-type quark masses and charged lepton 
masses m^ ., and is defined as 



I jM (MlMlm 2 w ) 



Ml 



77% - VTLo 



a=l,2 J 



oo 

dx I tdt 

o 



t + M? 



x In 



m^fl — x) + m^x + tx(l — x) 



mj.(l — x) + m^x + tx(l — x) 



In the asymptotic limit, i.e., in the limit of M 2 2 ^> m 2 



djti m w-> 



(15) 



the integral I jk i(Ml, M$,m%r) 



becomes independent of the flavor indices j, k, I and takes the form 



I jkl (Mt,Mlm 2 w ) 



Ml 



T 



Ml In 



M.j 



Ml In 



Mi 



1 [Ml In 



2 Ml 



Ml In 



2 Ml 



M{ 



Ml 



Ml - M, 



(16) 

When the index I is not equal to 3 (corresponding to no top quark inside the loop), this 
asymptotic expression well approximates the exact integral. The integral Ijkz corresponding to 
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top quark inside the loop obtained numerically is shown in Fig. 2 as a function of M 12 . Here 
we have taken the running top quark mass m t (fi = 1 TeV) = 150 GeV. It will turn out that 
will enter into the dominant source for neutrino masses. 

We now proceed to write down the neutrino mass matrix explicitly. First note that if the 
integral Ijki in Eq. (12) is flavor universal (as happens in the asymptotic limit as given in Eq. 
(16)), or if Ij k3 is the only dominant contribution, M v can be written as 

M v = rh I [YD d V T D u F*D e + D e F*D u VD d Y T ] . (17) 

where I is the flavor universal value of 1^. Now, this equation can be expanded in power series 
in light fermion masses. Combining constraints from flavor changing processes, we find, to an 
extremely good approximation, that the neutrino mass matrix is given by 

/ \^xy \y 

lm *xy ^xz \z+\^x I. (18) 



Here we have used the following definition of parameters 



2 m T & m T 2 2 m T 



X 



w 



TP* 
r 23 


, y 


= ^3 


z = 


Y23 


TP* 


^33 


Y33 


TP* 
r 32 


Y32 


( m c \ 


( m s s 




TP* 
r 33 


^33 


\m t ) 




Ijk3 



m = 2m F33F33 I jh3 . (19) 

As expected, the leading contributions are proportional to mtmbm T , see the definition of m in 
Eq. (14). One additional contribution, denoted as w in the (3,3) entry of Eq. (18), can be 
important, for a very restricted range of parameters. This is in spite of the additional mass 
suppression factors (m c / m t )(m s / m^) that appears in w, and occurs for the parameter choice 
-^33^33 < 10~ 5 and Mi )2 ~ (300 — 500) GeV. The Yukawa coupling combination F32I32 appearing 
in w is allowed to be order one, these couplings do not lead to excessive flavor violating processes 
even when the leptoquarks are light. (F 32 couples to cr, and Y 32 couples L 3 to sr, both 
of which are found to be rather un-constrained.) Contributions not shown in Eq. (18) have 
magnitudes 0.01 eV, because of mass suppression factors and the constraints on the Yukawa 
couplings arising from rare processes mediated by the leptoquarks. 

We shall consider the case w 1 for the most part of our analysis, but we shall also examine 
the special case where w ^> 1, which will be realized when the coupling F 33 takes anomalously 
small value and the leptoquark mass is less than about 500 GeV. 

The zero in the (1,1) entry of Eq. (18) is highly suppressed. We find, using the lowest 
non- vanishing contribution from Eq. (17) for this entry to be 

^ - » (f ) fe) '"»• (20 > 
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Now, unless (F* 3 /F 33 ) = (F* 3 /F£ 3 ) x takes values as large as m T /m e ~ 4000, this entry will be 
negligible for neutrino masses. That possibility is however inconsistent with an acceptable fit to 
neutrino data on the one hand, and flavor violating constraints on the other. To see this, let 
us note that a good fit to neutrino data, discussed in more detail later in this section, requires 
I ^33 1 ~ |^23| ~ |^i3| and x = I -F23 1 / 1 -F33 1 ~ iTir/nT-n ~ 16. Further, the leptoquark mass should 
not exceed about 10 TeV, or else the neutrino mass m 3 will be smaller than 0.05 eV, which 
would be inconsistent with the atmospheric oscillation data. For the same reason, F 33 <C 1, 
while allowed, would require the leptoquark mass to be much below a TeV. Then flavor violation 
constraints become important, fi — e conversion in nuclei sets an upper limit on the products 
I-F13-F23 1 < 10~ 4 and IY1.3Y23I < 10~ 4 corresponding to a leptoquark mass of 300 GeV. When 
these limits are combined with the need to obtain the right magnitude for m 3 , we see that |F 23 | 
cannot be smaller than 10 -2 . Consequently the (1,1) entry of M U) given in Eq. (20), is much 
smaller than 0.01 eV. This of course means that contributions to neutrinoless double beta decay 
from light neutrino exchange is negligible. See however, other contributions to this process via 
the exchange of a leptoquark and W boson in Fig. 10 discussed later, which turns out to be 
significant. 

The other remarkable feature of the mass matrix in Eq. (17) is that, although it may not 
be obvious from its form, M v has nearly zero determinant when w < 1. This would imply that 
one of the neutrinos is essentially massless in the model for most of the parameter space. This 
can be seen by observing that only the top quark mass, and not the c and u quark masses, has 
entered into M v when w is set to zero. In this case, each of the two terms in Eq. (17) has rank 
one, and the sum of the two terms has rank two, leaving the determinant to be zero. The first 
non-zero entry to the determinant in our expansion is found to be 



Again, with \x\ ~ m r /m M ~ 16, \y\ ~ 1, and w <C 1 (or equivalently 132-^3*2 lO 4 !^^), the 
determinant will be much less than (0.01 eV) 3 and thus m x ~ to a high degree of accuracy. 
Note that F 3 * 3 cannot be less than 10~ 3 , or else the neutrino mass scale will be too small. 

2.3 Predictions for neutrino oscillations when w <C 1 

We turn now to the predictions of the model when w < 1, which is realized is much of the 
parameter space. Since M u in Eq. (18) has only four (complex) parameters, there are two 
predictions for neutrino masses and mixings. These are summarized below: 



Det(M I/ ) 




(21) 



mi ~ 0, tan 2 #i 3 




12 



(22) 
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Figure 3: The absolute values of x, y, 



Here we have used the standard parametrization of the PMNS matrix. Furthermore, the two 
Majorana phases a and (3 (the phases of masses m x and m 2 respectively) are given by 

/3-25 + Ti, a~0, (23) 

with a ~ being a consequence of m,\ ~ 0. These predictions have been obtained in the context 
of textures before (see Ref. [13]). Here we have derived them without resorting to textures. 

To check the consistency of these predictions, we keep Am 2 olar = 7.65 x 1CT 5 eV 2 , Am 2 tm = 
2.40 x 1CT 3 eV 2 at their central values [14], and then vary sin 2 ^ = 0.304lo;o 22 ; in its allowed 
range. This gives from Eq. (22) 

sin 2 13 = {0.051,0.049, 0.046, 0.044} , (24) 

where the numbers correspond to the central value, and (1, 2, 3) a deviation in sin 2 9\2- This 
should be compared with the limit sin 2 ^ < (0.040, 0.056), at the (2, 3)<r level [14]. We see 
broad agreement, although sin 2 # 13 should be very close to the current limit. 

We can infer the values of (x, y, z) in Eq. (18) from neutrino data. For this purpose we use 
the central values from neutrino data: Am 2 olar = 7.65 x 10~ 5 eV 2 , Am atm = 2.40 x 10~ 3 eV 2 , 
along with sin 2 # 23 = 0.50, sin 2 # 12 = 0.304, and sin 2 13 = 0.051 [14]. We obtain 

12.21 - 4.17e i5 _ -0.112e i5 (2.927 + e i6 ) _ 0.473 + 0.342e i5 + 0.062e 2iS 

0.726 + 0.248e i5 ' V ~ -0.473 + 0.018e il5 + 0.062e 2i(S ' Z ~ 0.473 - 0.018e i<5 - 0.062e 2i<5 L 
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The absolute values of x, y, z are given in Fig. 3. We see that \x\ >> 1 and \y\, \z\ ~ 1. These 
values will directly affect the branching ratios of the leptoquarks, making a strong connection 
between collider physics and neutrino physics. 



2.4 A special case with w ^> 1 

It is interesting to see the constraints on the model when !o»l. This occurs for a small range 
of parameters when the leptoquark mass is less than about 500 GeV, and when I-F33Y33I <C 10~ 5 . 
The latter condition can be satisfied when either |F 33 | or \Y 35 \ is much smaller than 10~ 5 . Both 
of these couplings cannot be simultaneously very small, or else the (2,3) entry of M u will become 
unacceptably small. 

When F33 — > 0, we have the condition (M^)i 3 ~ 0, see Eq. (18). Since (M u )u ~ in this case 
we reproduce a well known texture model [16, 17]. This limit of our model leads to the following 
predictions for sin 2 6*13 and the phase parameters: 

t^w+s) * t^y(«-.fl) (27) 

cot^ #12 — tarn 9 12 cos 2 (a — p) 
This gives 0.012 < sin 2 #13 < 0.017, which is within reach of the next generation long baseline 
neutrino oscillation experiments. This case requires the leptoquarks to be lighter than about 500 
GeV, which are within reach of the LHC. 



2.5 The case of w ~ 1 



For the case when w ~ 1, there are no restrictions on the neutrino oscillation parameters, except 
the one arising from the vanishing of the (1,1) entry of M v . This relation can be taken to be a 
prediction for mi, and an additional relation for the phase parameters. These are given by 

m 3 



tan #13 ~ — sin #12 



1 + T± cot 2 12e ^-P) 



25 







IT 



cot~ 



m 2 

1 + a cos(/3 — a) 



(28) 



where a = — cot 2 9 



m,2 



12- 



a sin(/3 — a) 

In this scenario 9\ 3 can take any value between zero and its current 



experimental upper limit. For w 
scenario is testable at the LHC. 



1, the leptoquark mass cannot exceed about 500 GeV, so this 



2.6 Limit on the parameter //, 

There are certain restrictions on the parameter /i that enters in the neutrino mass Lagrangian of 
Eq. (6). It is this term that is responsible for the mixing of w -1 / 3 with x -1 ^ 3 . Since the SU(2) L 
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partner of w -1 / 3 does not mix with any other field, the \i term will induce new contributions to 
the electroweak p parameter (or equivalently the T parameter). This is because mixing via the 
\i term splits the masses of the SU{2) L doublet members. This mass-splitting must obey [12] 



^(AM) 2 < (57GeV) 2 , 



(29) 



where C = 3 for leptoquarks. In our model, apart from the \x term, there is a quartic scalar 
coupling, V D n\(Q^H)\ 2 , which contributes to the w -1 / 3 mass, and not to the co 2 ^ 3 mass, causing 
a further splitting. The parameter (AM) 2 in Eq. (29) is given by 

kv 2 + Ml - ml 



AM 



2Mi 



(30) 



where we have assumed that the two states are nearly degenerate. The two-loop induced neutrino 
mass is maximized when the mixing parameter sin 29 = 1. For this choice, one has the relations 
Mf 2 = ml =F fJ>v. When the leptoquark mass is well above 200 GeV, the contribution kv 2 in Eq. 
(30), which is at most of order (200 GeV) 2 , can be neglected in relation to the p:v term. In this 
case we obtain from Eq. (29) an upper limit on /i: (i < 0.65Mi. We shall use this constraint in 
deriving an upper limit on the leptoquark mass. 

There are other reasons why \i cannot be arbitrarily large: the positivity of leptoquark squared 
mass (needed for color conservation) and the perturbativity of the theory. The first of these 
conditions implies, from Eq. (10), that 



m,„m^ 



(31) 



The second constraint emerges because \i induces negative quartic couplings for the three fields 
X 1//3 , w 1 / 3 and H°. These induced couplings cannot exceed the corresponding tree-level quartic 
couplings which can at most be of order one. Otherwise the theory would be non-perturbative, 
or electric charge and color would break in the minimum of the theory. We parameterize the 
effective quartic couplings of the three fields as 

-C cS = A cff {H°*H ) 2 + A^ ff ( X *x) 2 + Kk (u*ou) 2 , (32) 

with H° = h/y/2. The 1-loop correction to these quartic couplings are shown in Fig. 4. 4 By 
evaluating these diagrams we obtains for the effective couplings, 



A c ff 

Ks 



32tt 2 
1 



m: 



m 2 J 2 

4 



m 2 + m 2 



m: 



m; 



^ln^ 

ml ml 



- 2 



128tt 2 [m 2 x - m 2 ) 2 



mi 



1287T 2 



mr. 



m 



2\2 
h) 



m 
m 



ln- 



m 



mi 



mi 



m 



* - 2 

2 Z 



In 



mi 



mi 



m 



h 



m 



h 



(33) 



4 For the h 4 term, there is a similar diagram generated by the SM (H^H) 2 quartic coupling. But here we focus 
on the diagrams generated by the fj, term. 
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Figure 4: The diagrams leading to the quartic couplings correction generated by \x. 



where the factor 3 in A e fr is a color factor. Since these effective couplings are negative, they must 
be smaller in magnitude compared to the corresponding tree level ones, otherwise the potential 
would become unbounded, suggesting instability of the vacuum. By demanding all these induced 
couplings to be less than one, so that the theory is perturbative with a stable vacuum, one can 
derive limits of /i. Of course, if the limit derived from Eq. (31) turns out to be more stringent, 
that will supersede the present limits. 

For illustration, we consider a few simplified cases. 

• m h pa <C m x . For m h = = 200 GeV and m x = 1 TeV, the bound from Eq. (31) 
gives a more stringent limit, i.e. /i < 1.1 TeV. For m x = 2 TeV and the other masses the 
same as above, the more stringent limit comes from the \ e s, viz., \i < 1.9 TeV. 

• m h <^ m w 7^ m x : In this case 



where m is the smallest of m w and m x . For example, if m u = 1 TeV, m x = 2 TeV, and 
m h = 200 GeV, /i < 5.7 TeV. 

• mh m^j pa m x : In this case 



where £ = max {1, ln^ — 1}. The choice of £ depends upon the value of m w and m x . For 
example, if m u = m x = 1 TeV and m h = 200 GeV, one uses £ = 1, so /i < 5 TeV. On 
the other hand, if m u = m y = 2 TeV and mh = 200 GeV, one uses £ = ( ln^ 2 - — 1 ) , so 
li < 9.4 TeV. 




(34) 



IM < m x (64vr 2 /£) 1/4 , 



(35) 
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2.7 Limit on leptoquark masses 

The scaling of the neutrino mass with the LQ mass goes as m v ~ {m t m b m T jj,v) / Mf . Since our 
neutrino mass matrix has a normal hierarchy structure, the neutrino mass scale mo should be 
around 0.02 eV. We shall assume that all the Yukawa couplings are bounded by 1, a requirement 
for perturbative treatment of the problem. As we have seen, a fit to the neutrino oscillation 
data requires F 33 /F 2 3 ~ m M /m r ~ 1/16. The maximum value of m is realized for sin 29 = 1, 
corresponding to largest allowed LQ mass. From these considerations we arrive at a limit on the 
LQ mass of 12 TeV. If the leptoquarks are in the multi-TeV range, all three members should be 
nearly degenerate, from the requirement of sin 29 ~ 1 and the p parameter constraint. A slightly 
better limit can be derived on their masses from the rare process \i — e conversion in nuclei (see 
the next section for details). From the experimental limit on this process, one can derive the 
upper limit on lY^l^l as a function of co 2 ^ 3 mass. Since neutrino mass fitting requires to be 
of the same order for % — 1 — 3, one can also determine an upper limit on Y33. Combining all of 
these, we obtain an upper limit of 10 TeV on M\. 

There are lower limits on leptoquark mass from the Tevatron where they could be pair pro- 
duced. The D0 and CDF experiments [15] have obtained limits on the second and the third 
generation leptoquarks of 316 GeV and 245 GeV respectively [15]. Since those experiments look 
for b-jets and missing energy, these bounds are strictly applicable for leptoquarks coupled to 
neutrinos. In our case, to~ 1 ^ 3 and x ^ 3 have direct couplings to the neutrinos (the latter from 
mixing) for which the quoted limits apply. As for w 2//3 , being a member of the same SU(2)l 
doublet as co^ 1 ^ 3 , it should be nearly degenerate with co" 1 ^ 3 , and thus the limits of Ref. [15] 
should apply to co 2 ^ 3 as well. 

3 Constraints from rare processes 

Since this model features lepton flavor (as well as total lepton number) violation, it is important 
to see how its parameters are constrained by the experimental data, especially those arising from 
rare decay processes that are forbidden in the SM. Here we derive a variety of limits on the 
couplings Yij and Fij as functions of the LQ masses. These limits should be (and have been) 
satisfied in the neutrino fits. The processes we shall consider are fi~ — > e~7, fi~ — > e + e~e~, 
fi — e conversion in nuclei, r~ — > e~r), t~ — > /i~r), B s d — B s d mixing, K — K mixing, D — D 
mixing, Df — > t^v decay, muon g — 2, n + — > fi + T^ decay, and neutrino less double beta decay. 
The expected rates for several of these processes are found to be in the interesting range for the 
next generation of experiments. 
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Figure 5: The diagrams leading to jj, — > cy. and S represent dj (uj) quarks and w 2 / 3 (x 1 ^ 3 ) respectively. 



3.1 — >■ e 7 

The rare decay — > e~7 arises from a one-loop diagram. In the limit of m e = 0, there are only 
three diagrams for this decay shown in Fig. 5. In general, the amplitude for this decay process 
can be written as [18] 

M = -e—u{ Ve )a^{f MX + /bi7 5 )«GvM9) (36) 

where m M is the muon mass, q = p M — p e is the photon momentum, and e M is the photon 
polarization vector. The effective couplings /mi and f E1 are found to be (repeated indices 
implies summation) 

,2 



3mi fY*Y 2j FuF* \ 

- ^(^3(,,)-^ 4 (^.)) 



3ml {Y*Y 2j . . / V,/-2 ; ,., ., 



where = m 2 d ./mj., x u . = m 2 ,/m 2 ; m u . and are the j-th generation of up- and down-type 
quark masses respectively; m w and m x are w 2//3 and x 1//3 masses respectively; and -F 3 (:r) and 
F±(x) are the dimensionless functions given by [19, 20] 

„ /N x (l-x)(5 + x) + 2(2 x + l) lux 



12 (l-x) z 
1 (1-i)(5i + 1) + 2x(2 + x) lni 
12 (1 -x) 



F 4 (x) = ^-^n^-r^-r^-r^ ^ (3g) 



These expressions assume that the leptoquark mixing angle 9 (see Eq. (9)) is zero. For non-zero 
9, one can make the following replacement: 



ml ^ MV 

X. a=l " 



where Ci = sin 2 9, ( 2 = cos 2 9, and M a are defined in Eq. (10). 
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Process BR 



Constraint 



. . i (3r\i / 1 O v 1 n-H I ^3(^)^23 1 2 | ^F 11 F* 1 + ^F 12 F| 2 +F 4 (x f )F 13 F 2 *3 | 2 3 , lxl0 -i9 

// ->■ e7 < 1.2 x 1U ^ 1 jjjj < GeV 4 



«r_K„v /11vin-7 I F «^) y !3 y 33 1 2 |A^ll^l + ^^32+^(x t )F 13 F*3| 2 16xlQ - 

r^e7 < i.i x ru ^4 I ^4 < GcV 4 

rr _v |F3fa)y 2 * 3 y 3 3| 2 | ^i^ + ^F^+F^QF^ | 2 6 , 7xl0 

r ->■ /i 7 < 4.5 x 10 jjjj + ^ < GcV 4 



Table 1: The constraints from ^ -> £jj. 



It is an excellent approximation to set the first and second generation quark masses to zero. 
The branching ratio for this decay is found to be 



27a f\F 3 (x b )Y{ 3 Y 23 \ 2 \^F U F* ± + ^F 12 F* 2 + F,(x t )F* 3 F 23 f 



167rG^ \ to* 

where a = e 2 /47r is the fine structure constant. Constraints from other ^ — > ^7 processes are 
presented in Table 1, where we have used the approximation = 0. 

An interesting feature of this analysis is that the couplings are only weakly constrained 
from these processes. This is owing to a GIM-like cancelation for the amplitude of this process 
from the two diagrams (Fig. 5). In the limit of down-type quark mass being zero, the two 
graphs exactly cancel. This cancelation occurs because the charge of the internal leptoquark 
(2/3) is twice as large and opposite in sign compared to the charge of of the internal down quark 
(—1/3). The amplitude of the graph where the photon is emitted from the scalar line is twice 
smaller compared to the one where the photon is emitted from the fermion line, leading to the 
cancelation. The non-canceling contribution is suppressed by a factor to^/to 2 in the amplitude, 
which weakens the limit. 

3.2 n~ — y e + e~e~ 

These processes also occur at the one-loop level in our model (see Fig. 6). Here, there exist 
photon, Z-boson, and the SM Higgs-boson exchange diagrams, as well as box diagrams. The SM 
Higgs diagrams are suppressed by the electron mass, and thus highly suppressed and therefore 
ignored. Therefore, only the photon, Z-boson, and the box diagrams will be evaluated. 

In the case of photon exchange, because the photon is now off-shell, and with the electron 
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Figure 6: The diagrams leading to /i — > 3e. qj and S represent dj (uf) quarks and w 2 / 3 (x 1 ^ 3 ) respectively. 



mass set to zero, the amplitude is given by 

e 2 

M 1 = —[u(p el )Y(fE0 + fM0l5)u(p^][u(p e2 )-f ll v(p e3 )} 



m. 



+ e 2 V — [u(Pel)^ U (fMl + /bi75)«(Pm)] K^h^O ^)] 

~ (Pel <-»• Pe 2 ) , (41) 

where g = p^—p e i- P^ and p e j are the incoming muon and outgoing electron momenta respectively. 
The couplings /eo and /mo are found to be 

3m 2 , (Y* 3 Y 2i FvF% 



ho ~ 32^ V^T* + "K"^' 1 ' 



where (? 1)2 = £(2 - 31n^), /i 1>2 = i(5 - 12m^), (? 3 = F^), and /i 3 = F 2 (x t ). The 
dimensionless functions -Fi(x) and F 2 (x) are given by [20] 

-4 + 9a; - 5a; 3 + 2(2x 3 + 3x - 2) lnx 
~ 36(1 -x) 4 ' 

(x - 1)(10 + sQe - 17)) + 2(s 3 + 6a; - 4) Ins 
F2(X) " 36(1 -x) 4 • (43) 

For the Z boson exchange diagram, the leading contribution after ignoring terms suppressed 
by m 2 /M 2 , with M being leptoquarks mass, is 

G 

M Z = -J= [VLUZ(Pel)'fu L (p lJl ) +T] R UR{Pel)l' X U R (p^\ 
X [u(Pe2h^(9v - 9Al5)v(Pe3)} 

~ (Pel ^ Pe2) , (44) 

where g v = — \ + sin 2 6w and g e A = — |. Since we set the first and second generation quark 
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masses to zero (which is a good approximation), the effective couplings, i] L and r] R , are given by 

Vl = ^Y* 3 Y 23 F 5 (x b ), 

m = ^F 13 F; 3 F 5 (x t ), (45) 

with 

. . x 1 — x + In x , , , 

= -2 (46) 

For the box diagrams, the amplitude can be written as 

•Mbox = ^J^{YY^ n (YYl) 21 [w(Pei)Yu L (p,)}[^ 

+ 64^ ^ F ^ 11 ( FFt ^ 2 [^M^M [^(Pe2)7^«(Pe3)] 

- (Pel <-» Pe 2 ) • (47) 

Here we have set all the quark masses to zero. 

We now proceed to the calculation of the total amplitude. In the loop integral functions, 
F1...F5, since the couplings are not constrained by yT — > e~7, one needs to include all the 
diagrams mediated by co 2 ^ 3 boson, since the Y^ couplings may be of order one. However, since 
the Z-boson exchange contribution is suppressed by the the bottom quark mass, it can safely 
be ignored. On the other hand, the Z exchange is significant for x 1 / 3 mediated process with the 
top quark inside the loop. The electron mass has to be included in the calculation in order to 
avoid infrared singularity. The branching ratio is found to be 



/ — + __, l 3V2 V f C^ k Y*Y 2] Y lk Y 2 l Cf k F* 3 F 2j F lk F* k 
BR(/i e+e e ) = — — - 3 — + -I °— 



\32tt 2 G f J \ mt 



(48) 
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where we have introduced Hermitian parameters C^ R , given by 



1 

7776 



7 2e 4 In 2 ^ - 108 (3e 4 + 2e 2 (YY*) J In ^ 
+ (449 + 68rr 2 ) e 4 + 486e 2 (FF t ) 11 + 243 [YY^)\] , for j, k = 1,2; 



s 288 



54e 4 F 1 (x 6 ) - 4e 2 (Ge 2 ^) + (1^%) (i n =| + itt) + 36e 2 {YY^) U Fi(x b ) 
+ 9e 2 (yy f ) n + 9 (yy^ij , for, j = 1,2 and A; = 3; 



32 



24e 4 F 2 (x fe ) + 8e 2 (YY^) U F 1 (x b ) + (YY^J , for j = k = 3; 



(49) 



7776 



( 288 In 2 ^| - 1584 In ^| + 108 In ^| + 272vr 2 + 2759 

V mj. rrig 

+ S ( FFt ) n ( n - 4 ln f|) + M (^) n > ^ j, k = 1, 2; 

[llF 2 (x t ) - 29F 4 (x t ) - 4F A (x t ) In g - 4 (F 2 (x t ) - 2F 4 (x t )) (in =| + itt) 
+^ [4v^(6sin 2 ^-l)G F m£F 5 (a; t ) (-41n^|-i47r + ll) 
Cf fc = <J + {FF^ u (l8F 2 (x t )-36F 4 (x t )-4ln^-t4n + ll) 

+M ( FFt )n [l6v^sin 2 w G F m 2 F 5 (z t ) + (FF t ) n ] , for j = 1,2 and fc = 3; 



f F|(x t ) - 3F 2 (x t )F 4 (x t ) + 2F 2 (* t ) (In ^ - ^ 



+T I^fo) - 2F 4 (x t )] [4V2(6sin 2 ^ - 1)G F m 2 x F 5 (x t ) + {FF^) U ] 
+2(1-4 sin 2 e w + 12 sin 4 9 W ) G 2 F m\Fl{x t ) 
{ +V2 (FFt) n sin 2 9 w G F m 2 x F 5 (x t ) + ^ (FF+)^ , for j = k = 3. 

For heavy quarks inside the loop, the factors C^ R are in the agreement with [20]. 

Regarding the neutrino oscillation data, one needs the limit on the coupling constants that 
appear in Eq. (19). From the branching ratio above, in the assumption that there is no acci- 
dental cancelation among the various contributions, so that one can omit terms like Y* 3 Y 2 3Yj k or 
FizF^Fjk with j, k — 1, 2, one gets, for a LQ mass of 1 TeV, 



Y 13 Y 23 \ < 7.6 x 10- 



|F 13 F 2 * 3 | < 1.8 x 10- 3 



if |F 13 | <C e and 



FisFo,! < 1.3 x 10" 



(50) 
(51) 



if IF 
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(d) 



Figure 7: The diagrams leading to fi — e conversion, qj and S represent dj (uj) quarks and cj 2 / 3 (x 1 ^ 3 ) 
respectively. 

3.3 11 — e Conversion in the nuclei 

Another lepton flavor violating process that may occur in this model is \x — e conversion in nuclei 
(fiN — y eN). The interesting type is the so-called coherent conversion in which the nucleus does 
not change its initial state [21]. In this section we shall discuss the implications of our model 
in such processes. It turns out that there exist tree level diagrams in addition to the one-loop 
diagrams (Fig. 7). These two sets of diagrams probe different Yukawa couplings. The loop 
diagrams involve the same couplings that appear in the neutrino mass matrix, and has to be 
given special attention. These loop diagrams are similar to fx~ —y e + e~e~ except the pair of 
up/down quark is attached to the photon and Z lines. Following [18], the branching ratio of this 
process is given by 

BR(M^eiV) = ^ Ee ^ff^ Fl \9l(A + Z) + g d L (2A - Z) + 2ZA, L | 2 

+L R, (52) 

where \p e \ and E e are the momentum and energy of the outgoing electron respectively, Z is the 
atomic number of the nucleus, A is the mass number of the nucleus, Z c g is the effective atomic 
number defined in [22], F p is the nuclear matrix element defined in [18, 21, 22], is the muon 
capture rate of nucleus N, and g£ d R are defined as 

9l = Vl Qu sin 2 9w^j , 9r = m Q - Qd sin 2 9 W 

n ( l O^B \ I ^ Y21 3y^ (yyt) ll( yyt) 21 

9r = VR t - Qu sin 6 W + 



v 4 ^ a vv J Am 2 x G F 32tt 2 G f Am 2 x 
Ag L = 2 ^^J {fEo(-ml) + f M i{-ml) + f M o(-ml) + /gi(-m 2 )) , 

AdR = 2 <^m 2 (f E °(~ m l"> + fMi{-rnl) - / M o(-m 2 ) - /gi(-m 2 )) . (53) 
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Element BR 



Constraint 



48 



Ti 



< 4.3 x 10 



-12 



17.61 0.53 



™2 I ™2 



+ 



< 



5.2x10' 
GeV 4 



208 



Pb < 4.6 x 10 



-li 



33.81 0.15 



' "T ^,2 



™2 ~r __2 



< 



9.7x10- 
GeV 4 



Table 2: The summary of fi — e conversion for Ti and Pb. The values of Z e g and Fp are taken from [21] whereas 
Tm are from [23]. 



Here Q u and Qd are the charges of the up- and down-quarks. As before, we can neglect terms 
proportional to t)l and F 3 (xd 3 ). The summary of constraints from this process is given in Table 
2, with 

af = (2A - Z) [327r 2 5 lj - {YY^ U ] , af = 2Ze%, bf = 2Ze 2 g v 
bf = (A + Z)[32n 2 5 1 ,-(FF^ u }+2Ze 2 F 4 (x Uj ) 

+ 8V2G F m 2 x F 5 (x t ) (\A- Z sin 2 9 W ) 5 3j , (54) 

where gj = gj\ q 2 = _ m 2^ and hj = hj\ q 2 = _ m 2. Here we take \p e \ ~ E e ~ m M . 

For 48 Ti, we get for LQ mass of 1 TeV, \Y* 3 Y 23 \ < 4.6 x 10~ 3 and |Fi 3 F* 3 | < 1.9 x 10~ 4 which 
are slightly stronger than the constraints from /i" — > e + e~e~. The tree-level diagrams alone give 
\Y*! Y 21 \ < 3.3 x 10" 7 and |F 11 F 2 * 1 | < 3.3 x 10" 7 . A new generation experiment called COMET 
has been proposed to reach a better sensitivity of 10~ 16 [24]. There is also discussion of testing 
fi — e conversion at a future Fermilab experiment. The prospects for these experiments observing 
new physics are good within our model. This is true even if the MEG experiment [25] obtains 
negative results for the /i — > cy decay. 5 

We can use the bound from /i — e conversion and neutrino mass fitting to predict the minimum 
branching ratio for fx — > 3e. Such a correlation is not possible for fj, — > cy because of the GIM-like 
cancelation that occurs there, and the fact that F\j is not constrained by neutrino data. Note 
that, Y i3 (assuming that only the bottom quark mediates the process) cannot be smaller than 
10~ 4 for LQ mass of 300 GeV, or else the induced neutrino mass will be too small. Therefore, 
the smallest branching ratio for \i — > 3e can be predicted which is presented in Fig. 8. We see 
that /i — > 3e can be substantial for a significant part of the parameter space. 

5 There is a class of model where [i — e conversion is log-enhanced compared to /i — > e-f, see Ref. [26]. 
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M ffl2/3 (TeV) 

Figure 8: The smallest branching ratio for \x — > 3e as suggested by neutrino mass fitting and ^ — e conversion. 

3.4 r~ — > e r/ and r~ — >• // r/ decays 

It turns out that some of the potentially large entries in the neutrino mass matrix are eliminated 
from the constraints from the decays r~ — > e~r] and t~ — > pTt). These processes can be mediated 
by tu 2 / 3 leptoquark at tree level. Here we ignore f} — rf mixing and work in the SU(3) flavor 
symmetric limit. The branching ratios for these processes are given by 



1 \Y 12 Y 32 \ 2 fim 3 T 



5127T 5 



mi 



r 



total 



1 \Y 22 Y 32 \ 2 fim 3 - 



5127T 5 



BR(r" e~rj) 

BR(r~ ->■ [i~r]) 

We use f v = 160 MeV and data from [12] to derive 

I^i2^3 2 | < 1-2 x KT 2 ( , . 

1 12 321 V300 GeV7 ' 



total 



(55) 



-V 



\Y 22 Y 32 \ < 1.0 x KT 2 ( 7 



(56) 



,300 GeV 

3.5 Muon anomalous magnetic moment 

The muon anomalous magnetic moment receives correction from the leptoquark exchange in our 
model and is given by 

2 



8{g-2) li = 3.9 x 10 



-10 



300 GeVj 





F 2 1 


2 


-F-22 


2 


F 2 3 




( 




+ 




- Ft(x t ) 


') 




1.0 




1.0 


1.0 





(57) 
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Again, for the same reason as in the case of // — > e^y, couplings are practically have no effects 
here. By comparing the new contributions with the 3<r anomaly in the experimental value, 
5(g — 2)^ xp = (24.6 ± 8.0) x 10~ 10 [27, 28], we see that this model can reduce or eliminate the 
discrepancy for order one values of the couplings F 2 i and F 22 . 



3.6 Upsilon Decay 

Lepton number violating T decay can set limit on the couplings Y. The decays we consider are 
T — > fi + iT and T — > fi + r~. The former is generated by the s-channel photon exchange and the 
latter is generated by the t-channel u 2 ^ 3 exchange. In our model the ratio of these two branching 
ratios is given by 

br(t^+t-) _ i \y* 3 y 33 \ 2 (m T \ 4 

BR(T^ n+fji-) 32 1 67r 2 a 2 \m w ) ' 1 ' 

where mj is the T mass. The upper limit of this ratio is 2.5 x 1CT 4 [29]. Therefore, we get less 
constrained value 

3.7 D~ — y t~v decay 

Another interesting feature of the leptoquarks is that they can mediate a lepton number violating 
D s decay: D~ — > l~v. Currently there is a discrepancy about 2a level between the experimental 
value and theoretical predictions based on lattice evaluation of the decay constant fo 3 - The LQ 
exchange can resolve the discrepancy. This issue has been studied in the framework of MSSM 
with broken i?-parity [27, 30] and in extensions of the standard model with leptoquarks [31]. It 
turns out that in this model this process occurs at tree level mediated by X a ^ 3 leptoquarks in 
addition to the SM process. Since the SM diagram leads to antineutrino final states, while the 
X a mediating processes result in neutrinos, there is no interference between the two. 6 We get 



BR(£>, = ^D s f D Gl\V ct 







2 


128MfG 2 F 


V cs 


2 




l+ .J.'TJH, ,. sin 2 2g(l-M 

x I 1-^1 <> (60) 

where rriD s and te> s are the D s mass and lifetime respectively while V cs is the (2,2) element 
of CKM matrix. Comparing results from lattice calculation, f^f = (241 ± 3) MeV [33], and 
the experimental result, = (257.5 ± 6.1) MeV [34], we see that at present there is a 2a 
discrepancy. If this discrepancy prevails, there is a natural explanation in the present context. 



6 Analogous lepton number violating muon decay has been explored as an explanation of the LSND anomaly 
in Ref. [321. 
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Assuming maximal mixing among the two charge —1/3 leptoquarks, we obtain the constraint 
(corresponding to the central values of the decay constants and f^ p ) 



Mi \V 450 GeV 



^^{rntev) (61) 

There is a similar contribution to the decay n~ — > fi~v e . From this process we find 

l^ifJil 2 . 2„„ {, M^ 2 



^'^"'"-^m^^^-wj ■ (62) 



By using experimental limits [12], one can derive 

i^<-(^) 2 (^)- 

3.8 New CP violation in B s — B s mixing 

The leptoquarks of the model mediate meson-antimseon mixing via one-loop box diagrams shown 
in Fig. 9. There is special interest in this process for the B s system, as there are experimental 
hints for non-zero CP violation in this system. We find that the present model can nicely explain 
the new CP violation that would be needed in B s — B s mixing, without generating unacceptable 
mixing in the B®, K° and D° systems. 

The new effective Lagrangian for B s — B s mixing in the model arising from the box diagrams 
of Fig. 9 is given by 

CT = {J * 1%R) {J * lM ' (64) 

The B s — B s transition amplitude is given by 

(B s \-CT\B s ) = m B JlB^(n)^(n), (65) 

where wb s and fs s are the B s meson mass and decay constant respectively, while B^ s (/j,) and 
r/f s (/i) are the bag parameter and the QCD correction factor evaluated at the scale /i ~ m b ~ 5 
GeV and their numerical values are 5f s (5 GeV) = 0.86 and r?f s (5 GeV) = 0.80 [35]. 

Recently the D0 collaboration has reported an excess in the like-sign dimuon asymmetry [36] 
defined as 

/y"++ _ AT 

A b sl = N++ - jpz = "(9.57 ± 2.51 ± 1.46) x 10" 3 , (66) 

where N ++ (N ) is the numbers of events containing two b hadrons that decay semileptonically 
into two positive (negative) muons. Eq. (66) can be written as a linear combination of two 
asymmetries [37, 36] 

A b sl = (0.506 ± 0.043)af, + (0.494 ± 0.043)a^, (67) 



23 



b co 2/s s,d b ej s ,d 



s,d "^~ /3 b s,d e) b 

Figure 9: The box diagrams of B s - B s Mixing. 
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a% = 1% Q ^J'S^ _Zi - (68) 



where a q sl (q = d, s) is defined as [36] 

T(B q -» /2 + X) - T(B q -» f£X) 
T(B q ^ l 2+X) + T(B q ^fi-X) 

In the SM, a d sl = -4.8t^xl0~ 4 and a s sl = (2.1±0.6) x lO" 5 [38], so that (A b sl ) SM = -2.3^ xl(T 4 
which is 3.2o" away from the current measurement (see Eq. (66)). A likely explanation is that 
there is a new source of CP violation in B s — B s mixing, which can arise from leptoquarks, as 
given in Eq. (65). 

The leptoquark contribution will modify the mass difference AM S in the B s system [39]: 

AM, = AMf + AM s new = AMf 1 1 + h s e 2ia ° \ . (69) 
In our model, the new contribution h s is given by 



**** = sUs m °-ww#v>- w 

The index i = 1 — 3 is to be summed in Eq. (70), but in practice not all three lepton flavors can 
simultaneously be significant, owing to lepton flavor violating constraints. Therefore we focus on 
the i = 3 contributions. 

A global fit to the CKM mixing parameters, including the new D0 data on dimuons, finds 
[39] {h s ~ 0.5, a s ~ 120°} or {h s ~ 1.8, a s ~ 100°} as the best fit values. The leptoquark 
contributions will be maximal when | V32 1 ~ 1, which is allowed in the model, and when the 
leptoquark mass is the smallest (m w ~ 300 GeV). In this case < 10~ 2 for % — 1, 2 in order to 
satisfy lepton flavor violating constraints. A fit to the neutrino oscillation data requires all the 
Y i3 (i — 1 — 3) to be of the same order with II33/Y23I ~ 3 corresponding to m w = 300 GeV. In 
this case, [i — e conversion in nuclei sets a limit II33I < 0.078. Using this upper limit in Eq. (70), 
and using |Mf 2 ^| = (9.0 ± 1.4) ps _1 [40], we find h s ~ 0.42. The phase a s is unconstrained, and 
can take the desired value of 120°. These are in excellent agreement with the global fit of Ref. 
[39]. 

We need to point out that the leptoquark model could affect the absorptive part of the mixing 
amplitude [41] : 7 

r 12s = r^ + r-r, (71) 



,2 



7 We thank Amol Dighe for useful discussion on this point. 
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where r°f™ is given by [41] 



\2 



rST = msJlB^^mlF^), (72) 

where -F(^) is a kinematic function, as a function of mj./ml. The numerical values of this 
function are F(£i) ~ 1 for £i = e, fi and -F(r) ~ 0.65. By using the same parameter values as for 
AM S , we see that the new contribution can modify the SM value of the width difference by at 
most 2%. There is currently a 2a discrepancy in Ar s , which should be explained as a statistical 
fluctuation within our model. 

Similar expressions hold for the mixing of other meson systems: 

(B d \-CT\B d ) = ggg mB d f 2 Bd B*<Mn?<M, 

(K\-CT\K) = ^^rn K f K B^M^), 

(D\-CT\D) = ^^^fnBiiMM- (73) 

The 7]i factors are found by following the procedures given in Ref. [35] to be 7/f d = 0.80, rff = 
0.76, ??f = 0.78, with Bf(5 GeV) = 0.86, B? (2 GeV) = 0.66, (2.8 GeV) = 0.865. 

The limits on those couplings appearing in K — K, D — D, and B^ — B^ mixings can be derived 
by using m K = 0.498 GeV, f K = 0.16 GeV, m D = 1.9 GeV, f D = 0.2 GeV, m Bd = 5.3 GeV, and 
f Bd = 0.24 GeV. They are found to be 

i y <^i < 96 x 10 " (ra) • 
i^i < 7 - 3 x 10 " (ra) ■ 

^ < L6 x 10 " (ao^v) ' < 74 > 

These constraints are all consistent with acceptable generation of neutrino masses and mixings 
within the model. 



3.9 Neutrinoless double beta decay 



The light neutrino exchange contribution to neutrinoless double beta decay vanishes in our model, 
owing to the zero of the (1,1) entry in M v . However, it can proceed via the vector-scalar exchange 
process, which does not need a helicity flip of the neutrino [42]. The relevant diagram is shown 
in Fig. (10). The effective Lagrangian of that process after Fierz rearrangement is 



-eflf 



e«7 M (l-7 5 )d 
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u{\ + 7 5 )rfe7^(l - 75)-7e c + -^ua a ^(l + 75)^7^(1 - ^)-a afi e c 



, (75) 
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Figure 10: The neutrinoless double beta decay diagram generated by this model. 



where q is the internal neutrino momentum. The first term is the scalar-pseudoscalar current 
and the second one is the tensor current, and the parameter e is defined as 



where 9 is the leptoquarks mixing angle. Following Pas et al. [43], the constraint on jV^Fnl 
is calculated by constructing nuclear matrix element in pn-QRPA model [44] and applying the 
result of OvftP half life, Ti/2(0z//3/3) > 1.2 x 10 25 yr, obtained from Heidelberg-Moscow experiment 
[45]. A straightforward calculation gives 



Here, both the scalar-pseudoscalar and tensor currents are taken into account. One sees that it 
is possible to observe neutrinoless double beta decay in the near future, even with hierarchical 
neutrino masses. 

4 Collider signals 

At hadron colliders, leptoquarks can be produced in pairs via quark-quark or gluon-gluon fusion 
and/or in association with lepton via quark-gluon fusion. The former depends on the QCD 
coupling constant and the latter depends on the \Fn\ 2 and l^iil 2 - The associated production may 
be important if the F and Y are as large as electromagnetic coupling e [46] . Since these couplings 
are constrained by rare processes such as neutrinoless double beta decay, it seems unlikely that 
single production will be dominant in our model. Therefore, here we focus on pair production of 
leptoquarks. 




(76) 




(77) 
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In w 2 / 3 pair production, each will decay to a charged lepton and a down quark with opposite 
charges. The resulting final state is two leptons and two jets. At the LHC with y^s = 14 TeV 
and mass of leptoquark equals to 500 GeV, the cross section is about 340 fb [46], decreasing to 
4.2 fb for a LQ mass of 1 TeV. The dominant background for this process is final states from 
Z + jets and tt production. For high leptoquark mass, only Z + jets is dominant. LHC can 
probe this leptoquark mass up to 1.2 TeV with integrated luminosity of 300 fb -1 . 

In Xa 1 ^ 3 pair production, there are two possible signals. First, each Xa 1 ^ 3 decays to up-quark 

— 1/3 

and a lepton, leading to two lepton and two jets signal. Secondly, one X a decays to charged 
lepton and up quark and the other one decays to neutrino and a down-type quark, leading to one 
lepton plus jets plus missing energy. The first signature is similar to the u 2 ^ 3 case. The second 
one has dominant background coming from W + jets and tt production. In Ref. [46], it has 
been shown that LHC can probe this leptoquark mass up to 1.2 TeV with integrated luminosity 
of 300 fb" 1 . 

The decays of the u?! 3 leptoquark into e + b, fi + b, r + b will occur in the following ratios: 

T{u 2 ' 3 -+ e + b) : T{u 2 ' 3 -+ fi + b) : T{u 2 / 3 -+ r + b) = \y\ 2 : \z\ 2 : 1 . (78) 

Measurement of any one of these branching ratios will determine the CP violating phase 5 via 
Eq. (25). This can of course be cross checked in long baseline neutrino oscillation experiments, 
especially since sin 2 #13 is large. If two branching ratios in the decay of co 2 ^ 3 are measured, that 
will determine the phase 5 in two different ways, allowing for another consistency check. 

The leptoquarks Xa 1 ^ 3 , (a = 1,2) which are linear combinations of x^ 1 ^ 3 an d w _1//3 will decay 
into charged leptons in the following ratios: 

rpr-vs fi-t) : r(x-V3 ^ T - t) = ^ . 1 (79) 

Note that this result holds for both X a ^ 3 , independent of the (co^ 1 ^ 3 — x^ 1 ^ 3 ) mixing angle 6. 
Measuring this branching ratio will determine providing another check for the model. Since 

— 1/3 

numerically \x\ ^> 1, we expect that at least one of the the X a will have [i~ in the final 

— 1/3 

state dominantly. It should be noted that the decay X a — > e~t will be suppressed, owing to 
constraints from fi — > cy. Of course, these Xa 1 ^ 3 fields also decay into Vidj. The lepton flavor 
composition in this mode would be near to impossible to determine. 

5 Conclusion 

In this paper we have presented a new two-loop neutrino mass generation model. This involves 
TeV scale leptoquarks, which can be directly tested at the LHC. The structure of the model is 
such that the neutrino mixing angle sin 2 #13 is predicted to be close to the current experimental 
limit of 0.05. The neutrino oscillation parameters are closely linked to the decay properties of 
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the leptoquarks. By measuring the branching ratios of the LQ bosons, in this model, one is 
measuring CP violating phase 5 of neutrino oscillations. These measurements also provide a 
number of cross checks by which the model can be falsified. 

We have discussed various rare decay processes mediated by the LQ bosons. The process 
— > e~7 may very well be suppressed, owing to a GIM-like mechanism, while the decay 

— > e + e~e~ and /x — e conversion in nuclei are within reach of the next generation experiments. 
A very interesting feature of the model is that the present hint of new physics contributions to 
CP violation in B s — B s mixing fits rather well here, without conflicting other well-measured 
meson mixings. Neutrinoless double beta decay can occur via LQ exchange, even when the 
neutrino spectrum has a normal hierarchy. Near future experiments will have a lot to say on the 
consistency of the model. 
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